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ABSTRACT 


As  the  U.S.  Navy  develops  new  technologies  which  enhance  automation  and 
reduce  crew  size  onboard  naval  vessels,  unmanned  vehicles  will  become  increasingly 
valuable  in  conducting  maritime  operations.  Effective  launch  and  recovery  systems 
(LARS)  are  necessary  for  unmanned  vehicles  to  efficiently  conduct  operations  at  sea. 
The  Towed  Body  is  a  LARS  with  a  wide  range  of  applications  for  unmanned  vehicle 
operations.  The  Towed  Body  can  be  evaluated  as  a  small  vessel  with  horizontal  and 
vertical  control  surfaces.  Since  it  is  being  towed,  the  directional  stability  of  the  Towed 
Body  requires  unique  consideration  due  to  the  presence  of  the  towing  force.  This  thesis 
examines  the  effect  of  varying  the  longitudinal  location  of  the  vertical  control  surfaces,  as 
well  as  the  effective  aspect  ratio,  size,  and  number  of  vertical  control  surfaces.  The 
results  identify  critical  stability  values  for  the  various  fin  configurations.. 
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I.  INTRODUCTION 


A.  SHIPBOARD  USE  OF  UNMANNED  VEHICLES 

As  the  U.S.  Navy  develops  new  technologies  which  enhance  automation  and 
reduce  crew  size  onboard  naval  vessels,  unmanned  vehicles  will  become  increasingly 
valuable  in  conducting  maritime  operations.  These  vehicles  will  be  capable  of 
performing  a  variety  of  missions  such  as  surface  search,  force  protection,  maritime 
interdiction,  and  mine  warfare. 

In  order  to  effectively  employ  unmanned  vehicles  from  onboard  ship,  efficient 
handling  systems  are  necessary  for  the  launch  and  recovery  of  these  vehicles.  The 
process  of  launching  and  recovering  unmanned  surface  vehicles  (USV)  and  unmanned 
underwater  vehicles  (UUV)  can  be  a  difficult  and  labor  intensive  evolution  under  adverse 
conditions.  High  winds  and  high  sea  states  produce  unfavorable  relative  motions 
between  the  unmanned  vehicle  and  ship,  increasing  the  risk  of  damage  and  personnel 
injury. 

B.  EXISTING  METHODS  FOR  LAUNCH  AND  RECOVERY 

The  process  of  launching  and  recovering  a  small  boat  from  a  ship  is  a  common 
procedure.  Current  naval  vessels  utilize  boat  davits  in  order  to  accomplish  this  task.  This 
system  requires  numerous  crew  members  to  operate  the  davit  and  hold  steadying  lines 
while  the  boat  is  raised  or  lowered  into  the  water.  Under  ideal  conditions,  this  is  a 
relatively  safe  and  easy  evolution.  High  winds  and  sea  states  cause  a  significant  increase 
in  the  relative  motion  between  the  ship  and  small  boat,  which  in  turn  increases  the  threat 
of  damage  to  the  vessels  and  for  personnel  injury. 

Some  smaller  U.S.  naval  vessels  and  U.S.  Coast  Guard  vessels  utilize  a  stem 
ramp  for  launch  and  recovery  of  small  boats.  This  method  of  launch  and  recovery  makes 
it  possible  to  conduct  the  evolution  with  fewer  personnel,  but  the  risk  of  damage  in  higher 
sea  states  still  exists  due  to  the  relative  motion  between  the  boat  and  the  ship.  Also,  the 
added  effect  of  operating  directly  in  the  wake  of  the  ship  must  be  considered  for  the 
performance  of  the  small  boat  attempting  to  maneuver  into  the  stern  ramp.  Overall,  as 
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compared  to  the  boat  davit,  the  stern  ramp  provides  a  safer  and  more  effective  method  for 
launch  and  recovery  of  small  boats.  [Sheinberg,  Cleary,  Beukema] 


Figure  1.  Stem  Ramp  on  USCGC  MATAGORDA 
[http://www.zodiacmilpro.com/product/sbytAI 

Originally  constructed  as  a  110  foot  Island  class  Patrol  Boat,  figure  (1)  shows  the 
USCGC  MATAGORDA  after  completing  a  service  life  extension  program  (SLEP).  This 
overhaul  lengthened  the  ship  to  123  feet  and  replaced  a  boat  handling  crane  with  a  stem 
ramp  [Gourley],  Several  U.S.  Navy  Cyclone  class  Patrol  Craft  received  similar  overhauls 
where  the  boat  handling  crane  was  removed  and  replaced  with  a  stern  ramp  for  small  boat 
launch  and  recovery. 


Figure  2.  USS  SQUALL  (PC  7)  with  boat  handling  crane  aft 


Cranes  are  utilized  by  some  naval  vessels  and  commercial  vessels  for  launching 

and  recovering  small  boats  and  underwater  vehicles.  Oceanographic  research  vessels  use 

cranes  or  A-Frames  to  handle  various  types  of  underwater  vehicles.  Existing  crane 

systems  require  personnel  to  operate  the  crane  and  handle  steadying  lines  to  the  vehicle. 
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High  winds  and  high  sea  states  again  cause  large  relative  motions  between  the  ship  and 
the  vehicle.  This  makes  the  process  of  connecting  the  vehicle  to  the  crane  cable  difficult 
to  accomplish  as  a  result  of  the  magnitude  of  the  ship  and  vehicle  motions. 


Figure  3.  AUV  A-Frame  LARS  by  Bluefin 
Robotics  [http  ://www.  bluefim-obotics.com/launch.htm] 


C.  FUTURE  LAUNCH  AND  RECOVERY  SYSTEMS 

The  Innovation  Center  at  Naval  Surface  Warfare  Center,  Carderock  Division 
published  a  report  in  2003  assessing  new  concepts  for  unmanned  vehicle  launch  and 
recovery.  This  study,  titled  SCOUVO  (Surface  Combatant  Optimized  for  Unmanned 
Vehicle  Operations)  evaluated  three  different  hull  foims  along  with  several  launch, 
recovery,  and  onboard  handling  systems  optimized  for  use  with  unmanned  vehicles.  For 
launch  and  recovery,  the  study  examined  five  systems:  an  automated  crane  system 
(Homing  Crane),  Variable  Geometry  Cradle,  Towed  Body,  Paravane,  and  Chinese 
Lantern. 

The  Homing  Crane  allows  for  reduced  manning  when  operating  the  system  and 
provides  a  soft  initial  connection  between  the  ship  and  unmanned  vehicle.  The  crane  is 
motion  controlled  and  maneuvers  itself  to  pick  up  the  vehicle  alongside  the  ship.  This 
can  be  used  for  both  USV  and  UUV  operations. 
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The  Variable  Cradle  is  utilized  with  a  stem  ramp.  The  cradle  changes  the 
geometry  of  the  landing  area  to  provide  proper  support  for  the  particular  UV  which  it  will 
hold.  This  allows  the  cradle  to  hold  a  variety  of  vehicles  with  varying  hull  forms.  This 
system  is  also  USV  and  UUV  capable. 


Figure  4.  Variable  Cradle  Concept  [SCOUVO] 

The  Towed  Body  is  a  maneuverable  body  which  is  towed  astern  of  the  launch  and 
recovery  ship.  It  allows  for  a  soft  connection  between  the  ship  and  UV  and  can  be 
complemented  with  a  crane  or  stern  ramp.  This  system  will  operate  with  both  US  Vs  and 
UUVs. 

The  Paravane  is  also  towed  behind  the  ship  and  snares  the  UV.  This  concept  is 
based  on  the  A/N37U-1  Mine  Clearing  Set. 

Only  for  use  with  UUVs,  the  Chinese  Lantern  can  retrieve  multiple  vehicles  at 
once.  This  system  drags  a  line  behind  the  ship  with  several  hook-up  locations  along  the 
length  of  the  line.  The  UUVs  maneuver  to  the  hook-up  points  on  the  line,  and  the  ship 
retrieves  the  vehicles  when  desired  [SCOUVO]. 
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Figure  5.  Chinese  Lantern  UUV  Retrieval  [SCOUVO] 


D.  TOWED  BODY 

This  thesis  will  examine  the  directional  stability  of  the  Towed  Body  launch  and 
recovery  system.  Since  this  concept  has  the  potential  to  function  with  a  wide  variety  of 
vehicles  and  may  be  used  with  a  stem  ramp  or  crane,  it  is  a  possible  launch  and  recovery 
solution  with  many  applications. 

The  concept  for  the  Towed  Body  was  initially  developed  late  in  the  Cold  War  for 
underway  refueling  of  surface  combatants.  The  body  would  be  located  at  the  end  of  a 
refueling  line  and  allow  for  a  low  signature  underway  refueling  evolution  by  eliminating 
the  need  for  standard  topside  refueling  rigging  [Mulhern], 

The  Towed  Body  concept  was  reintroduced  as  a  UV  launch  and  recovery  system 
in  the  SCOUVO  report.  In  addition  to  enabling  a  soft  connection  between  the  UV  and 
ship  for  retrieval,  the  towline  can  also  carry  with  it  a  means  to  refuel  the  vehicle  or 
provide  a  connection  for  data  to  be  passed  between  the  UV  and  ship.  This  allows  for 
multiple  sorties  of  the  UV  without  needing  to  recover  the  vehicle  onboard  ship  each  time. 


5 


Figure  6.  Towed  Body  Concept  [SCOUVO] 


The  conceptual  design  of  the  towed  body  has  a  wingspan  of  approximately  three 
feet  and  a  body  length  of  about  five  feet.  The  control  surfaces  (wings  and  fins)  are 
moved  by  electrically  driven  actuators.  The  fins  are  located  on  the  bottom  of  the  body  to 
allow  the  hook-up  device  for  grasping  the  UV  to  be  located  on  top  of  the  body. 

The  retrieval  capacity  of  the  Towed  Body  system  depends  primarily  on  the  size  of 
the  tow  cable  and  the  capacity  of  the  constant  tension  winch  used  to  reel  in  the  tow  cable. 
The  system  is  envisioned  to  accommodate  any  UV  up  to  20,000  pounds,  approximately 
the  size  of  a  U.S.  Navy  11  meter  RHIB. 

To  accomplish  the  hook  up  between  the  Towed  Body  and  the  UV,  the  Towed 
Body  will  be  equipped  with  a  guidance  system  similar  to  that  of  a  heat  seeking  missile. 
The  UV  will  carry  an  optical  emitter  for  the  towed  body  to  locate  and  then  maneuver  to 
that  location.  This  is  an  experimental  guidance  system  that  has  not  been  tested  in  an 
operational  environment  [Mulhern],  The  harness  on  the  Towed  Body  used  to  connect  to 
the  UV  operates  like  a  mousetrap.  Once  the  UV  trips  the  harness,  the  Towed  Body  is 
reeled  in  back  to  the  ship.  The  vertical  fins  on  the  towed  body  are  able  to  retract  in  order 
to  allow  the  body  and  the  UV  to  be  pulled  onto  a  stem  ramp  without  damage. 

The  Towed  Body  design  differs  slightly  for  USV  and  UUV  use,  as  seen  in  Figure 
6.  This  study  will  focus  only  on  USV  operations.  The  Towed  Body  is  assumed  to  be 
operating  in  the  vicinity  of  the  ocean  surface  and  will  only  maneuver  in  the  x-y  plane. 
No  change  in  depth  will  occur.  This  assumption  permits  the  wings  to  be  ignored  and  as 
such  will  not  be  included  in  the  model  developed  in  this  thesis. 
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II.  EQUATIONS  OF  MOTION 


A.  EQUATIONS  OF  MOTION 

A  marine  vehicle  has  six  degrees  of  freedom  in  which  it  can  move  when  operating 
without  any  restraints.  Figure  7  shows  the  three  translational  degrees  of  freedom  (surge, 
sway,  and  heave)  and  the  three  rotational  degrees  of  freedom  (pitch,  roll,  and  yaw)  on  and 
about  the  x,  y,  and  z  axes. 


Figure  7.  Coordinate  Frames  on  USS  SQUAFF  (PC  7)  [after  USS  SQUAFF 

command  photograph] 

The  equations  of  motion  for  all  six  degrees  freedom,  as  derived  in  [Healey  (2001)]  are: 
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SURGE :  m -  vr  +  wq  - xG  (q2  +  r2 }  +  yG  (pq -  r )  +  zG  (/>r  +  =  X 
SWAY  :  rn^y  +  ur -  wp  +  xG  (pq  +  r)— yG  (p2  +  r2  )  +  zG  ( qr - p)J  =  Y 
HEAVE :  m[w-uq  +  vp  +  xG(pr-q)+yG(qr  +  p)  —  zG  (/r  +  ^2)]  =  Z 
ROLL:  IJ+^l-I^qr+I^pr-qyi^-r^-I^pq  +  r) 

+  m\_yG  ( w—uq+vp)—zG  (v+ur-wp) J  =  K  (2.1) 
PITCH:  lyq+(lx-lz)pr-Ixy(qr  +  p)+Iyz(pq-r)+Ixz(p2-r2) 

-  m\ycG  (w—uq+vp)— zG  {ii-vr  +  wqj^  =  M 
YAW:  Ij  +  ily-I^pq-I^-q^-I^pr  +  ^  +  I^qr-p) 

+  m  (v  +  ur  -  wp )  -  yG  (ti  -  vr  +  w’<y)]  =  N 

B.  YAW  AND  SWAY  SIMPLIFICATION 

The  analysis  in  this  thesis  is  focused  on  the  directional  stability  of  the  towed  body 
for  USV  operations,  so  several  simplifications  can  be  made  to  reduce  the  equations  of 
motion.  First,  any  environmental  forces  that  may  act  on  the  towed  body,  such  as  wind 
and  waves,  will  be  ignored.  For  USV  retrieval,  the  towed  body  will  operate  in  the 
vicinity  of  the  ocean  surface  and  will  not  perform  any  depth  changes.  Also,  it  is  assumed 
that  the  towed  body  will  be  subjected  to  a  constant  towing  force  and  that  it  will  have  the 
same  forward  velocity  as  the  towing  vessel.  These  assumptions  allow  the  equations  of 
motion  to  be  reduced  down  to  the  two  equations  for  sway  and  yaw.  Surge,  heave,  pitch, 
and  roll  are  neglected  in  this  study. 

Following  the  justification  and  simplified  derivative  notation  of  the  Society  of 
Naval  Architects  and  Marine  Engineers  (SNAME)  [PNA],  the  equations  of  motion  are 
reduced  as  follows: 


SWAY  : -Yvv +(m-Yi)v-(Yr-mU)r-Y.r  =  0 
YAW  :  -Nvv  -  N,v  -  Nrr  +  (/z  -  A, )  r  =  0 


(2.2) 


The  terms  Yvv  and  Yrr  represent  the  lateral  force  due  to  side  slip  velocity  and  yaw  rate 
respectively.  The  terms  Nvv  and  Nrr  represent  the  moment  caused  by  side  slip  velocity 
and  yaw  rate.  The  term  m  represents  the  mass  of  the  vessel,  U  is  the  forward  velocity, 
and  /  is  the  mass  moment  of  inertia  about  the  z  axis. 
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These  equations  are  nondimensionalized  in  order  to  permit  easier  use  of  available 
data.  Nondimen sional  terms  are  denoted  with  a  prime  superscript.  The  sway  equation 
consists  of  force  terms  and  is  nondimensionalized  by  dividing  with  ( p  1 2)  IsV 2 .  The 
yaw  equation  consists  of  moment  terms  and  is  nondimensionalized  by  dividing  with 
(/?/2)LV2.  This  yields: 


SWAY  :  -Y'V  +  {m' - Y'.  ) v' - (r/  -  m'jr'  - Y'r’  =  0 
YAW  :  -Nv'v'  -  N-v'  -  N'r'  +  (//  -  N-  )  r1 '  =  0 


(2.3) 


Forward  velocity,  U,  is  nondimensionalized  with  ship’s  speed,  so  u'  =1  and  the  term  is 
dropped  from  the  equation. 

Additional  sway  forces  and  yaw  moments  are  produced  by  the  deflection  of  a 
vertical  control  surface,  such  as  a  rudder.  Figure  8  shows  the  turning  moments  produced 
by  rudder  deflection  on  a  ship. 


Figure  8.  Rudder-induced  turning  moments  [PNA] 


The  rudder  deflection  angle  is  identified  by  SR  ,  and  Ys  and  Ns  are  the  derivatives  of  Y 

and  N  with  respect  to  the  rudder  deflection  angle.  These  terms  may  also  be 
nondimensionalized,  and  are  included  in  the  sway  and  yaw  expressions  as  follows: 
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(2.4) 


SWAY  :  -Y'V  +  {rn'  - Y.’ J v'  -  (y  V)  r'  -  Y'r'  =  Y'S5R 
YAW  :  -Nv'v' -  N-'v'  - N'r' + -  N  '  jr'  =  N'sSR 

For  the  purpose  of  this  study,  the  forces  and  moments  produced  by  the  rudder  will 
be  considered,  but  only  with  the  rudder  held  at  zero  degrees. 

C.  CONTROLS  FIXED  STABILITY 

When  a  vessel  is  moving  in  a  straight  line  path  and  encounters  an  instantaneous 
disturbance,  the  subsequent  path  of  the  vessel  illustrates  its  motion  stability.  If  the  vessel 
resumes  a  straight  line  path,  regardless  of  direction  after  the  disturbance  occurs,  it  is 
considered  to  have  straight  line  stability.  A  vessel  with  directional  stability  is  able  to 
resume  the  straight  line  path  in  the  same  direction  as  the  original  path.  Figure  9 
illustrates  the  different  types  of  motion  stability. 


ORIGINAL  STRAIGHT 


\  straight  line  stability 


* 


ORIGINAL  straight 
—  uN£  Path 

n  DIRECTIONAL  stability 

_  (WITH  eOfctfLEJl  StAAlUfY  IfvOEXESj 
ORIGINAL  STRAIGHT  & 

line  path 

HI  DIRECTIONAL  STABILITY 

[WITH  P£AL  STAJS5LTY  INDEXES} 


original  strai5ht-$3* 


FINAL  PATH  IS  STRAIGHT  BIJT 
DIRECTION  CHANGED 


FINAL  PATHNAME 
DIRECTION  AS 


ORIGINAL  PATH  SvlT 
DIFFERENT  POSITION 


-  FINAL  RAT*,  SAME  AS  CASE  X 


FINAL  PATH,  SAME 
DIRECTION  AND  POSITION  AS 


ORIGINAL  FATH 


IY  POSITIONAL  MOTION  STABILITY 

INDICATES  INSTANTANEOUS  DISTURBANCE 

I  I  I  I  1 _ I _ _ 


Figure  9.  Motion  stability  [PNA] 

The  straight  line  stability  of  a  vessel  can  be  deteimined  with  the  characteristic 
equation: 


Act2  +  Bcr  +  C  =  0 


(2.5) 


where 
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(2.6) 


A  =  (I[-  N' )(«'  -  Y!)  -  (Y!  -  m’x’G  ){N'~  m'x’G ) 

B  =  -(rz-N'f)Y:-(m'-Y!)(N'r-m'x'Gu') 

-  (y; -  mV)  (N’  -  m'xG )  -  (Y! -  m'x'G )  N[ 
C  =  (N’r -m'x'Gu')Y;-{Y'r-m'u')N'v 


In  order  to  meet  the  stability  criteria,  A,  B,  and  C  must  be  greater  than  zero.  The  value  of 
C  is  termed  the  stability  index.  Generally,  a  positive  stability  index  indicates  that  the  ship 
is  stable.  Alternatively,  the  stability  index  can  be  found  in  normalized  maimer  with  the 
expression: 


n'(w;-m'4) 


(2.7) 


Equation  2.7  allows  for  easier  comparison  and  plotting  of  the  stability  index  [Papoulias], 
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III.  DEVELOPING  HYDRODYNAMIC  COEFFICIENTS 


A.  TOWED  BODY  COMPONENTS 

The  Towed  Body  LARS  is  currently  a  conceptual  design,  so  there  are  no  existing 
models  or  specifications  available  for  evaluating  the  directional  stability  of  the  design.  In 
order  to  develop  a  model  for  the  purpose  of  this  thesis,  existing  data  for  similar  body 
shapes  were  utilized  in  a  nondimensional  format. 

1.  Body  of  Revolution 

The  geometry  of  the  towed  body  can  be  modeled  with  an  existing  body  of 
revolution,  the  Defense  Advanced  Research  Projects  Agency  (DARPA)  SUBOFF  model, 
developed  at  the  Naval  Surface  Warfare  Center,  Carderock  Division  (NSWCCD).  This 
model  was  utilized  by  [Wolkerstorfer]  in  his  work  with  SLICE  pods.  Wolkerstoifer 
approximated  the  SUBOFF  body  as  a  body  of  revolution  with  an  elliptical  nose, 
cylindrical  mid  body,  and  a  conical  base  section. 


The  nondimensional  hydrodynamic  coefficients  for  the  SUBOFF  body  were 
determined  through  captive  model  tests  at  NSWCCD.  These  values,  shown  in  Table  1 
are  used  in  this  study  of  the  Towed  Body. 
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y: 

-0.005948 

y; 

-0.013278 

Yr 

0.001811 

y; 

0.000060 

K 

-0.012795 

K 

0.000202 

K 

-0.001597 

N' 

-0.000676 

Table  1.  SUBOFF  Hydrodynamic  Coefficients 


Nondimensional  values  for  mass  and  mass  moment  of  inertia  were  not  available 
for  the  approximate  SUBOFF  body.  For  an  acceptable  approximation  of  these  values,  the 
coefficients  from  the  U.S.  Navy  DSRV  were  used  to  establish  a  linear  relation  between 
mass  and  Y'  and  between  /'  and  N. .  This  relation  was  then  applied  to  the  Y.  and  N'- 

values  for  the  SUBOFF  body  to  obtain  acceptable  values  for  m'  and  /' .  \DSRV\ 


iri 

0.013594 

/: 

0.00084997 

Table  2.  SUBOFF  values  interpolated  from  DSRV  data 
2.  Vertical  Control  Surfaces 

The  fins  for  the  Towed  Body  are  modeled  using  the  standard  NACA  0015 
geometry  as  an  all-movable  rudder.  This  rudder  design  has  a  square  tip  shape,  taper  ratio 
of  0.45,  sweep  angle  of  zero,  and  a  thickness  to  chord  ratio  of  0.15.  This  standard 
geometry  scales  easily  with  various  aspect  ratios  and  chord  lengths. 


Figure  11.  NACA  0015  Rudder  Cross  Section  [PNA] 
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Figure  12.  Specifications  and  Terminology  for  an  all-movable  rudder  [PNA] 

The  geometric  aspect  ratio  is  determined  by  dividing  the  mean  span  by  the  mean 
chord  length.  The  effective  aspect  ratio  takes  into  consideration  the  groundboard  effect 
of  the  fin  being  mounted  to  a  flat  surface.  For  the  towed  body,  a  groundboard  condition 
is  assumed,  so  the  effective  ratio  is  two  times  the  geometric  aspect  ratio  [PNA], 

The  longitudinal  location  of  the  rudder  is  described  with  the  variable  xF ,  the 
longitudinal  distance  from  the  center  of  the  body  of  revolution  to  the  center  of  the  rudder. 
Positive  values  of  xF  indicate  that  the  rudder  is  located  forward  of  amidships,  or  the 
center  of  the  Towed  Body. 

Experimental  data  provides  the  most  accurate  method  of  determining 
hydrodynamic  coefficients.  Since  experimentation  is  not  always  practical,  theoretical 
calculation  is  acceptable.  The  theoretical  calculations  in  this  study  determine  the 
hydrodynamic  coefficients  of  the  fixed  vertical  control  fins  within  a  linear  range. 
Following  the  theoretical  development  in  Principles  of  Naval  Architecture,  Volume  III, 
the  hydrodynamic  coefficients  for  a  fixed  fin  are  expressed: 
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Af  =  area  of  fin 
b  =  geometric  span  (3.2) 


Using  Jone’s  Formula  for  low  aspect  ratio  fins,  =  ja ,  where  a  is  the 

effective  aspect  ratio  of  the  fin.  Since  the  groundboard  condition  is  assumed  for  the 
towed  body  design,  a  =  2aG.  These  relations  reduce  and  (7/)  from  equation  3.1 

to: 


B.  TOTAL  DERIVATIVES  FOR  BODY  AND  FINS 

In  order  to  develop  the  hydrodynamic  coefficients  for  the  towed  body  including 
the  fins,  the  individual  derivatives  for  the  body  and  each  appendage  are  summed.  For  a 
towed  body  design  using  only  one  vertical  fin,  the  hydrodynamic  coefficients  would  be: 
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J'MnO.+OV), 

K=(K)t+W), 

K'=(K')t+(K% 

K=(K)t+(Kl 

v =(n'),+ 07), 

*;=(*.vw), 

rz-OtVOV), 

w;=(w;),+(jv;), 

When  two  fins  are  evaluated  in  this  study,  both  fins  are  assumed  to  be  at  the  same 
longitudinal  location  on  the  body,  and  equally  offset  from  centerline.  From  equation  3.4, 
this  yields: 

r:=(r:)t+2(r:)f 

N',=(N:\+2(N:)f 

K=(N'r\+  W), 

JV=(iV)»+2(r;)/ 

N',=(N',\+2(K), 

r;=(Y;i+2(r;)f 

"l=W)>+2(W), 
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IV.  TOWING 


A.  GEOMETRY  AND  NOMENCLATURE 

The  towed  body  must  be  evaluated  in  a  slightly  different  manner  than  a  free 
maneuvering  ship  due  to  the  presence  of  the  tow  cable  force.  It  is  assumed  that  the 
hydrodynamic  coefficients  are  determined  in  the  same  fashion  as  previously  described, 
but  the  stability  criterion  takes  into  account  the  impact  of  the  towing  force. 

Figure  13  illustrates  the  geometry  and  terms  for  developing  the  equations  of 
motion  for  a  towed  vessel. 

u0  =  towing  vessel's  forward  velocity 
l  =  length  of  the  unstrained  towline 
T  =  tension  of  the  towline 
L  =  length  of  the  towed  body 

x  =  distance  from  towed  body's  CG  to  tow  connection 
ya  =  towed  body's  offset  from  towing  ship's  line  of  motion 
i//  =  towed  body  yaw  angle 
x  =  distance  between  two  vessels  in  x  direction 


Figure  13.  Towing  Geometry  [Bemitsas] 


The  addition  of  the  tow  cable  force  alters  the  relevant  nondimensional  equations  of 
motion  in  equation  2.3  as  follows  [Latorre]: 
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SWAY:  (Y!-  m')  %  +  Yv%  -jy0  +  Yfr'  -  (Y'vu'  -  T/)  ij,' 


f 

(  x'  ) 

Y’u'+T’ 

1+^ 

l  n 

I 

+ 

V 

j 

YAW  :  NJ0  +  NJ0  -  j  x'p  y’  +  (N{  - I'z )  y  -  {N\u'  -  N’r )  \j/’ 


(4.1) 


f 

r  y ) 

\ 

Ny+T'x'p 

i+-£- 

l  n 

+  Nlu' 

V 

J 

y/’  =  0 


B.  STABILITY  CRITERIA  FOR  TOWING 

The  characteristic  equation  used  for  evaluating  the  stability  of  a  towed  vessel,  as 
derived  in  [Latorre]  and  [Bemitsas]  is  of  the  form: 


Acr  +Bcr  +Ccr  +Dcr  +  E  =  0 


(4.2) 


Where 


A=(Y!-ni)(N';-i:yN'X 

B=(Y!-rri)(K -N\u')+y;{N',  -I'z)+(Y-u'-Y^)N',  -YJK 


T . 


c=Y:(K-W)—(KO+(r-nt) 


(  -Ol 

-N'u'-T'Y 

V  p 

l  n\ 

y;n'-Y;)N[+Y^Tx’^ 


V  r  y 


D=r 


-N'u'-TV 
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7Y 


+j(w -KHW- -ty-f+K 


f  -Ol 

Y'u'+T' 

V 

l  n\ 

T 

f  Y) 

(  -Ol 

E=— 

A/V+7Y 

— 

Y'u'+T' 

r 

v  p 

l  rJ. 

V 

l  n\ 

r 


(4.3) 


Towing  is  determined  to  be  stable  if  all  real  parts  of  the  roots  of  the  characteristic 
equation  are  negative.  The  work  of  [Latorre]  and  [Bemitsas]  further  show  that  B,  D,  E, 
and  BCD- AD2  -BE  must  all  be  greater  than  zero.  Also,  since  T ,  / ,  and  xp  are  the 

most  important  parameters  in  towing,  the  terms  in  A,  B,  C,  D,  and  E  may  be  organized 
with  subscript  ‘0’  for  teims  in  the  maneuvering  stability  criteria,  ‘1’  for  coefficients  of  T, 
and  ‘2’  for  coefficients  of  Tjl .  This  results  in: 
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A=A> 


B  =  B, 


o 


c=c0+r 


c 

Q+— 

1  v 


Co  =  Y'N'r  +  N'v  ( m' -  Y') 

C\  ={m'-Y-)xrp  +N'. 

C2  =-(N'f -rs)+(m'-Y!)X,p2+x,p  (Nl  +Y!) 
ZU 


D=T' 


D,+- 


V 


d{=n:-y;x'p 


d2 = -YX +N[x'p+ (y; -y:)x'p-  N'r + N't 


T' 

E=—E2 
E2  =  Dx 


(4.4) 


Ultimately,  there  are  two  active  criteria  for  ship  towing  stability.  The  first 
criteria,  Ri,  is  determined  by  the  following  inequality: 


N' 

R,  :  x'„  > 

1  p  yr 


(4.5) 


The  second  criteria,  R2,  is  determined  by  the  inequality 


R, :  T'>- 


(4.6) 


where 


ai  -  (fi0Q A  ~  \Dl ) + +  B0C2Dl  -  2AqDvD2  ) + —{boC2D2  -  AqD2  ) 

2  (4.7) 

a2  =  B(iC0D1  +—(boC0D2  -bj~e2) 

The  inequality  in  R2  can  also  be  used  to  identify  the  critical  value  for  T'  which 
makes  towing  stable  because  the  values  ax  and  a2  are  independent  of  speed,  resistance, 
and  tow  cable  tension. 


21 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


22 


Y.  COMPUTER  MODELING  AND  PROBLEM 
CONFIGURATION 


A.  STABILITY  EVALUATION 

The  calculations  conducted  in  this  thesis  will  determine  if  the  Ri  stability  criteria 
is  satisfied  for  a  range  of  fin  sizes,  effective  aspect  ratios,  and  longitudinal  locations  on 
the  towed  body.  Also,  the  critical  value  for  T' ,  the  nondimensional  tow  cable  tension, 
will  be  identified  for  the  various  fin  configurations. 

B.  VARIABLES 

The  four  variables  used  to  find  the  stability  characteristics  for  a  wide  range  of  fin 
configurations  are:  Number  of  fins  (n),  effective  aspect  ratio  (a),  nondimensional 
longitudinal  location  of  the  fin  ( x'F ),  and  the  nondimensional  root  chord  length  of  the  fin 

(CL). 

1.  Number  of  Fins  (n) 

Calculations  are  conducted  for  a  Towed  Body  with  no  vertical  control  surfaces, 
one,  and  two  vertical  control  surfaces.  When  one  fin  is  used,  it  is  located  on  the 
centerline  of  the  body.  For  the  case  of  two  vertical  control  surfaces,  the  fins  are  located 
at  the  same  longitudinal  location  on  the  body,  and  are  placed  symmetrically  about  the 
centerline  of  the  body. 

2.  Effective  Aspect  Ratio  (a) 

The  effective  aspect  ratio  is  two  times  the  geometric  aspect  ratio  following  the 
groundboard  assumption  previously  discussed.  Calculations  are  conducted  for  effective 
aspect  ratios  ranging  from  0.5  to  6.0  at  an  interval  of  0.5. 

3.  Longitudinal  Location  of  Fin  ( x'F  ) 

The  nondimesnional  value  x'F  is  found  by  dividing  xF  by  the  vessel  length.  The 
fins  will  be  evaluated  at  locations  ranging  from  amidships  to  the  stern.  As  a  rudder  is 
moved  forward  on  a  ship,  it  becomes  less  effective.  A  Towed  Body  configuration  with 
fins  forward  of  amidships  is  not  sensible  due  to  this  and  the  presence  of  the  towing  force 
near  the  bow.  Therefore,  the  values  for  xF  will  range  from  -0.04  to  -0.5  at  an  interval  of 
-0.02. 
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4.  Root  Chord  Length  (C'root ) 

The  value  for  the  root  chord  length  determines  the  overall  dimensions  of  the  fin 
because  the  aspect  ratio  has  already  been  assigned.  The  root  chord  length  is  also 
nondimensionalized  with  the  length  of  the  towed  body.  This  variable  will  range  from 
0.02  to  0.3  at  an  interval  of  0.01. 

C.  MATLAB  PROGRAMMING 

MATLAB  version  7.0.1  was  used  to  conduct  the  calculations  for  this  thesis.  In  an 
effort  to  make  the  program  easily  tailored  for  a  variety  of  variables,  it  was  developed  as  a 
number  of  small  programs  for  each  individual  component  and  type  of  calculation.  This 
simplifies  the  effort  required  for  the  user  to  change  the  specifications  for  the  body  of 
revolution  and  the  fins.  It  also  allows  additional  components,  such  as  the  wings,  to  be 
easily  added  to  the  program.  The  hydrodynamic  coefficients  are  calculated  in  an 
individual  file,  and  the  stability  criterion  are  also  calculated  in  an  individual  file.  A 
separate  file  is  used  to  run  the  entire  program  and  assign  the  range  of  values  for  each  of 
the  variables  through  the  use  of  several  ‘for’  loops.  The  nested  ‘for’  loops  produce  data 
for  a  total  of  16,705  rudder  configurations.  All  of  this  data  is  then  stored  in  an  array,  and 
additional  program  files  are  used  to  develop  the  plots  seen  in  Figures  15  through  37.  The 
m.files  developed  to  conduct  these  calculations  may  be  reviewed  in  appendix  A. 
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VI.  RESULTS 


A.  Rj  STABILITY  CRITERIA 

The  first  condition  evaluated  is  the  Ri  stability  criteria.  This  condition  is 
independent  of  the  tow  cable  tension,  Towed  Body  resistance,  and  towing  vessel’s  speed. 
The  data  shows  that  for  each  effective  aspect  ratio  evaluated,  there  is  a  critical  value  for 
Croat  a|U  x'f  f°r  which  the  Ri  criteria  is  satisfied.  Evaluated  with  no  fins,  the  Towed 

Body  does  not  meet  the  Ri  stability  criteria.  Table  3  summarizes  the  critical  values  for 
the  Ri  stability  criteria  for  the  fin  configurations  examined  in  this  study.  The  value 
identified  as  x'P  indicates  the  most  forward  position  on  the  body  which  that  fin  may 

r  CRITICAL  x  j  j 

be  located  and  still  meet  the  Ri  stability  criteria. 


a 

n=l 

n=2 

c 

root 

Xp 

p CRITICAL 

C 

root 

Xp 

r CRITICAL 

0.5 

<0.13 

All  fail  R, 

<09 

All  fail  Ri 

0.13 

-0.44 

0.09 

-0.48 

0.14 

-0.32 

0.10 

-0.30 

0.15 

-0.22 

0.11 

-0.18 

0.16 

-0.14 

0.12 

-0.06 

0.17 

-0.04- 

>0.12 

All  satisfy  Ri 

>0.17 

All  satisfy  Ri 

1.0 

<0.07 

All  fail  Ri 

<0.05 

All  fail  Ri 

0.07 

-0.32 

0.05 

-0.30 

0.08 

-0.14 

0.06 

-0.06 

>0.08 

All  satisfy  Ri 

>0.06 

All  satisfy  Ri 
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<0.05 


All  fail  R1 


<0.03 


All  fail  Ri 


1.5 

0.05 

-0.22 

0.03 

-0.48 

>0.05 

All  satisfy  Ri 

0.04 

-0.06 

>0.04 

All  satisfy  Ri 

<0.04 

All  fail  Ri 

<0.03 

All  fail  Ri 

2.0 

0.04 

-0.14 

0.03 

-0.06 

>0.04 

All  satisfy  Ri 

>0.03 

All  satisfy  Ri 

<0.03 

All  fail  Ri 

0.02 

-0.30 

2.5 

0.03 

-0.22 

>0.02 

All  satisfy  Ri 

>0.03 

All  satisfy  Ri 

<0.03 

All  fail  Ri 

0.02 

-0.06 

3.0 

0.03 

-0.04 

>0.02 

All  satisfy  Ri 

>0.03 

All  satisfy  Ri 

3.5 

0.02 

-0.32 

All 

Satisfy  Ri 

>0.02 

All  satisfy  Ri 

4.0 

0.02 

-0.12 

All 

Satisfy  Ri 

>0.02 

Satisfy  Ri 

>4.0 

All 

Satisfy  Ri 

All 

Satisfy  Ri 

Table  3.  Critical  Values  for  Ri  Stability  Criteria 


The  data  shows  that  as  the  effective  aspect  ratio  increases,  the  size  of  the  fin  can 
decrease  while  meeting  the  stability  requirement.  Each  step  increase  in  effective  aspect 
ratio  yields  a  smaller  fin  which  meets  the  requirement.  With  an  aspect  ratio  greater  than 
4.0,  all  of  the  modeled  fin  configurations  satisfy  the  Ri  stability  criteria. 
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The  data  also  demonstrates  that  increasing  from  one  to  two  vertical  control 
surfaces  has  an  impact  on  the  satisfaction  of  the  R\  criteria.  For  a  given  effective  aspect 
ratio,  the  two  fin  configuration  will  meet  the  Ri  criteria  with  a  smaller  fin  compared  to 
the  one  fin  configuration. 

B.  G'  STABILITY  COMPARISON 

The  Ri  condition  is  compared  to  the  G[  stability  index  for  the  same  fin 
configurations.  Although  the  G’h  equation  is  for  use  with  a  free  running  vessel,  it  has 
been  computed  for  the  Towed  Body  and  compared  with  the  Ri  results.  Table  4  shows  the 
G[  and  Ri  comparison  for  a=l  and  C’root  =0.08.  Due  to  the  vast  amount  of  data  compiled, 

this  small  range  of  data  is  shown  because  it  exhibits  the  typical  response  seen  over  the 
entire  range  of  data  where  the  transition  from  unstable  to  stable  occurs  for  both  Ri  and 
G[ .  As  found  by  evaluating  the  Ri  criteria,  a  towed  body  with  no  vertical  control 

surfaces  does  not  meet  the  G'h  stability  criteria. 


Number  of 

effective 

x'F 

C'root 

51(i) 

Ri  (2) 

RI  satisfied? 

Gh 

fins 

aspect  ratio 

(0=no,  1=yes) 

1 

1 

-0.04 

0.08 

0.48 

0.54316 

0 

-2.7233 

1 

1 

-0.06 

0.08 

0.48 

0.52848 

0 

-2.4462 

1 

1 

-0.08 

0.08 

0.48 

0.51381 

0 

-2.1611 

1 

1 

-0.1 

0.08 

0.48 

0.49913 

0 

-1.8749 

1 

1 

-0.12 

0.08 

0.48 

0.48445 

!  o 

-1.5938 

1 

1 

-0.14 

0.08 

0.48 

0.46978 

1 

-1.3231 

1 

1 

-0.16 

0.08 

0.48 

0.4551 

1 

-1.0667 

1 

1 

-0.18 

0.08 

0.48 

0.44042 

1 

-0.82736 

1 

1 

-0.2 

0.08 

0.48 

0.42574 

1 

-0.60675 

1 

1 

-0.22 

0.08 

0.48 

0.41107 

1 

-0.40557 

1 

|  1 

-0.24 

0.08 

0.48 

0.39639 

1 

-0.22377 

1 

1 

-0.26 

0.08 

0.48 

0.38171 

1 

-0.06075 

1 

|  1 

-0.28 

0.08 

0.48 

0.36704 

1 

0.084473 

1 

1 

-0.3 

0.08 

0.48 

0.35236 

1 

0.21311 

1 

1 

-0.32 

0.08 

0.48 

0.33768 

1 

0.3265 

1 

1 

-0.34 

0.08 

0.48 

0.323 

1 

0.42603 

1 

1 

-0.36 

0.08 

0.48 

0.30833 

i 

0.51308 

1 

i  i 

-0.38 

0.08 

0.48 

0.29365 

1 

0.58893 

1 

i 

-0.4 

0.08 

0.48 

0.27897 

1 

0.65483 

1 

i 

-0.42 

0.08 

0.48 

0.2643 

1  i 

0.71189 

1 

i 

-0.44 

0.08 

0.48 

0.24962 

1 

0.76113 

1 

i 

-0.46 

0.08 

0.48 

0.23494 

i 

0.80348 

1 

i 

-0.48 

0.08 

0.48 

0.22026 

i 

0.83976 

1 

i 

-0.5 

0.08 

0.48 

0.20559 

i 

0.87072 

Table  4.  Ri  and  G'h  Comparison  for  a=l,  C'root=  0.08 
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In  Table  4,  a  red  zero  for  Ri  indicates  that  the  Towed  Body  is  unstable  for  that 
given  fin  configuration,  hi  the  G'h  column,  a  negative  value  is  red  and  also  indicates  an 

unstable  configuration.  A  cell  is  highlighted  in  yellow  where  the  Ri  and  G'h  stability 
conditions  do  not  match. 

The  two  stability  criterion  generally  do  not  concur  with  which  configuration  is  the 
critical  value,  marking  the  change  from  an  unstable  to  stable  Towed  Body.  However, 
both  Ri  and  G'h  exhibit  similar  behavior  across  the  range  of  data,  where  the  Ri  criteria  is 
satisfied  with  a  fin  which  is  smaller  or  further  forward  on  the  Towed  Body  as  compared 
to  the  G'h  value.  This  trend  allows  the  G'h  values  to  be  used  as  a  means  to  validate  the  Ri 
findings. 

Figures  15  through  19  exhibit  the  effect  of  fin  size  and  location  on  the  value  of 
G[  for  the  range  of  effective  aspect  ratios  evaluated  in  this  study.  Each  figure  displays  a 
separate  plot  for  one  fin  and  two  fins  in  order  to  demonstrate  the  effect  that  this  variable 
has  on  the  value  of  G'h .  Points  on  the  curves  which  fall  between  zero  and  one  in  the 
figures  represent  a  stable  configuration. 

These  figures  show  that  as  the  root  chord  length  increases,  the  fin  can  be  placed 
further  forward  on  the  body  while  meeting  the  G[  stability  requirement.  Also,  as  the 

effective  aspect  ratio  is  increased,  smaller  sized  fins  satisfy  the  G'h  stability  requirement. 

In  addition,  increasing  from  one  to  two  fins  allows  smaller  sized  fins  to  meet  the  stability 
requirement,  and  for  the  same  sized  fins  to  meet  the  criteria  further  forward  on  the  Towed 
Body. 
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c.  r2  stability  criteria 

In  order  for  the  Towed  Body  to  be  evaluated  as  stable,  both  the  Ri  and  R2 
criterion  must  be  satisfied.  The  R2  stability  criteria  is  dependent  on  the  tow  cable  tension, 
which  is  a  function  of  the  Towed  Body  resistance  and  the  speed  of  the  towing  vessel. 
Since  this  study  is  using  a  conceptual  design  for  the  Towed  Body,  actual  tow  cable 
tension  values  are  difficult  to  accurately  calculate.  This  is  due  to  the  lack  of  adequate 
design  data  to  estimate  the  resistance  of  the  Towed  Body  at  various  speeds.  By  using 
nondimensional  values,  the  R2  criteria  can  be  evaluated  to  find  the  critical 
nondimensional  tow  cable  tension  ( Tcrit )  which  is  required  to  satisfy  the  stability  criteria. 

The  T'rlt  values  obtained  in  this  study  can  be  converted  into  a  force  if  a  Towed  Body 

model  is  developed  in  detail.  The  force  in  the  tow  cable  can  then  be  used  to  determine 
minimum  towing  speeds  required  to  meet  the  R2  stability  criteria. 

This  section  presents  numerous  graphs  showing  how  T'rit  is  affected  by  varying 
the  number  of  fins,  the  effective  aspect  ratio,  the  root  chord  length,  and  the  longitudinal 
location  of  the  fin.  In  any  case  where  a  value  for  T'ril  falls  at  or  below  zero,  it  can  be 

assumed  that  particular  Towed  Body  configuration  will  meet  the  R2  criteria  without  any 
tension  present  in  the  tow  cable.  It  is  not  possible  for  a  negative  tension  to  exist  in  the 
tow  cable,  as  the  physical  cable  will  not  permit  a  compression  force  of  any  significant 
magnitude  to  exist  within  the  cable. 

1.  Impact  of  n  on  T'rit 

Figure  21  demonstrates  the  effect  that  the  number  of  vertical  control  surfaces  has 
on  T'rit  for  effective  aspect  ratios  of  one  through  four,  with  C'root  held  constant.  For  a 

Towed  Body  without  any  fins,  T'rit  is  found  to  be  less  than  zero.  However,  this 

configuration  cannot  be  evaluated  as  being  stable  because  with  no  fins,  the  Towed  Body 
did  not  meet  the  Ri  criteria. 

Comparing  the  curves  for  one  and  two  fins,  T'rit  is  greater  for  the  two  fin 
configuration  at  any  given  longitudinal  fin  location.  The  graphs  in  Figure  21  also  show 
that  for  both  one  and  two  fin  configurations,  the  value  for  T'rit  decreases  as  the  fin  moves 
aft. 
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,C' root  =  0.16  a=2,  C' root  =  0.16 


Figure  21.  Effect  of  n  on  T'ril  for  a  =  1  to  4 
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2.  Impact  of  a  on  T'rit 
a-  Ku  vs.  x'F 

Figures  22  through  25  demonstrate  the  effect  that  the  effective  aspect  ratio 
of  a  vertical  control  surface  has  on  T'rit .  The  figures  show  the  effect  for  both  one  and  two 
fin  configurations. 

Analysis  of  the  graphs  in  Figures  22  and  23  shows  that  as  the  effective 
aspect  ratio  is  increased,  T'rit  will  also  increase.  Also,  as  the  fin  is  moved  aft  from  an  4 

position  of  approximately  -0.10,  the  value  of  T'rit  decreases  for  any  given  effective  aspect 
ratio. 

Comparing  the  graphs  in  figure  22  for  one  fin  and  figure  23  for  two  fins,  it 
is  also  apparent  that  a  higher  value  of  T'rU  is  calculated  for  two  fins  of  a  given  size,  aspect 
ratio,  and  location. 

b.  Tb  vs.  C'root 

Figures  24  and  25  also  show  the  effect  of  increasing  the  effective  aspect 
ratio,  but  in  this  case  T'rit  is  plotted  against  C'root  for  specified  values  of  4  •  Here  it  is 

seen  that  for  C'root  values  generally  less  than  0.07,  T'rit  falls  below  zero,  indicating  that 
the  R2  condition  is  satisfied  with  no  tension  in  the  tow  cable. 

As  found  in  the  previous  section,  an  increase  in  the  root  chord  length 
results  in  an  increase  in  the  corresponding  values  of  T'rit .  However,  the  graphs  in  Figures 

24  and  25  show  this  relation  only  holds  true  for  C'root  values  ranging  from  0.05  to  0.12. 
As  C'oot  continues  to  increase  beyond  0.12,  there  is  minimal  effect  on  T'rit  for  any  given 
effective  aspect  ratio. 

As  before,  an  increase  in  the  number  of  fins  produces  a  higher  value  of 
T'rit  with  the  same  root  chord  length  and  location. 
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Figure  22. 


Effect  of  a  on  Tfcrit  for  specified  C[oot  values,  n  =  1 
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Figure  23.  Effect  of  a  on  Vcrit  for  specified  Croot  values,  n  =  2 
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Figure  24.  Effect  of  a  on  T'rit  for  specified  x'F  values,  n  =  1 
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Figure  25.  Effect  of  a  on  T'rit  for  specified  x'F  values,  n  =  2 
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3.  Impact  of  C'root  on  T'rit 

The  impact  which  the  size  of  the  rudder  has  on  the  value  of  T'rit  is  shown  in 
Figures  26  through  31.  In  these  figures,  T'rit  vs.  x'F  is  plotted  for  specified  values  of 
C'root  and  a  for  both  one  and  two  fin  configurations. 

It  is  apparent  by  examining  these  figures  that  T'rit  decreases  as  the  fin  is  moved 
aft  from  an  approximate  x'F  position  of  -0.1.  As  x'F  approaches  zero,  or  the  center  of  the 
body,  the  T'rit  values  experience  a  considerable  shift  for  mid-range  values  of  C'root . 

In  general,  T'rit  increases  as  the  effective  aspect  ratio  increases.  Increasing  from 
one  to  two  fins  also  increases  the  corresponding  T'rit  values.  As  expected,  this  concurs 
with  the  results  from  previously  discussed  figures. 

Figures  26  through  31  exhibit  an  interesting  trend  for  T'rit  based  on  C'toot  and  a. 
In  Figure  27,  the  graph  for  a=1.5,  n=l  shows  that  the  upper  limit  for  T'rit  corresponds 
with  C'root  =  0.30.  As  C'root  decreases  from  0.30  through  the  first  three  step  intervals,  the 
T'rU  values  are  near  the  T'rit  values  for  C'rgot  =  0.30.  A  significant  shift  in  T'rit  values 
does  not  occur  until  C'root  decreases  to  0.14. 

Reviewing  Figures  28  through  31  for  n=l  and  a  increasing  to  6.0,  it  is  seen  that 
the  upper  boundary  of  T'rit  values  increases.  Also,  as  a  increases,  the  T'rit  values  for 

decreasing  values  of  C'wot  approach  the  upper  boundary  of  T'rit .  In  Figure  28  for  a=3  and 
n=l,  a  significant  drop  in  T'rit  does  not  occur  until  C'root  =  0.10.  The  graph  for  a=6  and 
n=l  in  Figure  31  shows  that  the  T'rit  values  for  C'root> 0.06  all  fall  along  the  upper 
boundary  for  T'rit . 

The  graphs  for  n=2  in  Figures  28  through  31  show  that  the  same  trend  occurs  for 
a  Towed  Body  with  two  fins.  The  difference  from  a  Towed  Body  with  one  fin  is  that  the 
upper  boundary  values  for  T'rit  are  greater.  Also,  the  T'rit  values  for  smaller  sized  fins 

approach  the  upper  boundary  of  T'rit  values  with  lower  effective  aspect  ratios  compared 

to  the  corresponding  configuration  where  n=l. 
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Effect  of  C'root  on  T'nt  for  a  =  0.5  and  1.0 
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Figure  27.  Effect  of  on  Tfcrit  for  a  =  1.5  and  2.0 
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Figure  28.  Effect  of  C'root  on  T'rit  for  a  =  2.5  and  3.0 
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Figure  29.  Effect  of  C'root  on  T'nt  for  a  =  3.5  and  4.0 
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Figure  30. 


Effect  of  C'root  on  T'rit  for  a  =  4.5  and  5.0 
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Figure  31.  Effect  of  C'root  on  T'rit  for  a  =  5.5  and  6.0 
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4. 


Impact  of  x'F  on  T'rit 
The  effect  of  moving  the  vertical  control  surface  longitudinally  on  the  Towed 
Body  is  presented  in  Figures  32  through  37.  Here,  T'rit  vs.  C'root  is  graphed  for  specified 

values  of  x'F  and  a.  The  figures  show  this  effect  for  one  and  two  fin  Towed  Body 
configurations.  The  graphs  in  these  figures  show  similar  trends  as  seen  in  previous 
figures. 

In  most  cases  where  C'rool  is  less  than  0.10,  T'ril  is  less  than  or  equal  to  zero, 

indicating  that  no  tow  cable  tension  is  necessary  to  satisfy  the  Ri  criteria.  Before 
assuming  stability,  the  fin  configuration  must  be  checked  against  the  Ri  criteria  from 
Table  3.  Where  C'rool  is  less  than  0.10  for  lower  effective  aspect  ratios,  the  Ri  criteria  is 
often  not  satisfied. 

As  demonstrated  in  previous  figures,  T'rit  increases  as  the  effective  aspect  ratio 
increases  for  C'root  values  greater  than  0.10.  Also,  T'nt  values  for  n=2  are  higher  than 
corresponding  fin  configurations  for  n=l. 

As  C'root  increases  above  0.15  up  to  0.30,  theTT  values  become  somewhat 
constant.  In  this  range  for  C'root ,  the  graphs  show  that  T'rit  decreases  as  x'F  increases. 
This  trend  corresponds  with  the  results  seen  in  Figures  26  through  31. 
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Figure  32. 


Effect  of  x'F  on  T'rit  for  a  =  0.5  and  1.0 
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Figure  33. 


Effect  ofx^  on  T'ril  for  a  =  1.5  and  2.0 
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Figure  34.  Effect  ofx^  on  T'rit  for  a  =  2.5  and  3.0 
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Figure  35. 


Effect  of  x'p  on  T'rit  for  a  =  3.5  and  4.0 
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Figure  36. 


Effect  of  x'p  on  T'rit  for  a  =  4.5  and  5.0 
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Figure  37.  Effect  ofx^  on  T'rit  for  a  =  5.5  and  6.0 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  study  evaluated  the  motion  stability  of  the  Towed  Body  LARS  with  a  variety 
of  fin  configurations.  To  evaluate  towing  stability,  the  Ri  and  Ri  stability  criterion  were 
assessed  for  each  fin  configuration.  The  Ri  criteria  was  also  compared  with  the  G'h 
stability  index  for  a  free  running  vessel. 

Examination  of  the  critical  values  for  Ri  in  Table  3  and  the  graphs  in  Figures  15 
through  20  finds  predictable  trends  in  meeting  this  stability  criteria. 

•  The  further  aft  a  fin  is  placed  on  the  Towed  Body,  the  more  likely  it  is  to  be 
stable. 

•  Increasing  fin  size  allows  the  fin  to  be  placed  further  forward  and/or  have  a 
lower  effective  aspect  ratio  while  meeting  the  stability  criteria. 

•  Increasing  the  effective  aspect  ratio  allows  the  fin  to  be  placed  further  forward 
and/or  be  smaller  in  size  while  meeting  the  stability  criteria. 

•  Adding  additional  fins  reduces  the  minimum  effective  aspect  ratio  and  fin  size 
necessary  to  meet  the  stability  criteria.  It  also  allows  the  fins  to  be  placed 
further  forward. 

Evaluation  of  the  R2  stability  criteria  yields  critical  values  for  the  nondimensional 
tow  cable  tension.  The  data  displayed  in  Figures  21  through  37  also  show  predictable 
trends  for  the  minimum  value  of  T  required  to  meet  this  stability  criteria.  The  following 
conclusions  are  for  fin  configurations  placed  between  the  stem  and  a  distance  of  10%  of 
the  Towed  Body’s  length  aft  of  the  center  of  the  body  (-0.50  <  x'F  <  -0.10 ). 

•  The  use  of  multiple  fins  of  the  same  size  and  effective  aspect  ratio  yields  a 
higher  T’crit  . 

•  An  increase  in  effective  aspect  ratio  will  increase  T'rit . 
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•  Increasing  the  root  chord  length  of  a  fin  while  holding  the  effect  aspect  ratio 
constant  has  little  effect  on  T'rU  for  a  root  chord  length  greater  than 

approximately  15%  of  the  Towed  Body’s  length  ( C'root  >0.15). 

•  T'rit  decreases  as  a  fin  is  moved  aft  from  a  position  of  x'F  =-0.10  towards  the 
stem. 

In  order  for  the  Towed  Body  to  possess  motion  stability,  both  the  Ri  and  R2 
criterion  must  be  satisfied.  While  a  smaller  size  and  lower  effective  aspect  ratio  fin 
design  permits  a  lower  nondimensional  tow  cable  tension  to  pass  the  R2  criteria,  it  is  less 
likely  to  satisfy  the  Ri  criteria.  However,  as  the  fin  is  placed  further  aft  on  the  Towed 
Body,  the  Ri  criteria  is  more  likely  to  be  satisfied  and  a  lower  nondimensional  tow  cable 
tension  is  necessary  to  satisfy  the  R2  criteria. 

B.  RECOMMENDATIONS 

There  is  a  great  deal  of  additional  study  to  be  conducted  on  the  design  of  the 
Towed  Body.  The  work  conducted  in  this  paper  only  addressed  motion  in  the  x-y  plane, 
and  focused  on  just  two  degrees  of  freedom:  sway  and  yaw.  Since  the  Towed  Body 
concept  is  for  both  USV  and  UUV  operations,  the  wings  should  be  modeled  and  motion 
stability  analysis  should  be  conducted  for  pitch  and  heave.  The  effect  of  environmental 
forces  on  the  motion  stability  may  also  be  addressed  in  future  work. 

Once  a  detailed  model  of  the  Towed  Body  is  developed  where  weights  and 
dimensions  are  known,  the  R2  criteria  may  be  used  to  determine  specific  operating 
requirements,  such  as  towing  speed  and  the  necessary  capacity  of  the  tow  cable. 
Subsequent  work  may  also  address  the  performance  of  the  movable  control  surfaces  for 
sway,  yaw,  roll,  pitch,  and  heave  motions.  This  will  in  turn  permit  study  of  the  guidance 
system  to  be  used  on  the  Towed  Body. 
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APPENDIX 


MATLAB  PROGRAM  FILES 


A.  CALCULATING  HYDRODYNAMIC  DERIVATIVES 

1.  Body  of  Revolution  Hydrodynamic  Derivatives  (bodyspec.m) 

%LT  Scott  D.  Roberts 

%Body  Specifications  for  Towed  Body  USV  retrieval  system 

%Reference  Principles  of  Naval  Architecture  Volume  III  and  SLICE  Thesis 
%20  January  2005 

%Body  Dimensions 

L  =  5;  %length  in  feet 

T  =  1.111;  %draft  in  feet  (actually  diameter  of  body  of 

rotation) 

xG  =  0;  %longitudinal  center  of  gravity 

%Hydrodynamic  Derivatives  for  Body  —  Nondimensional 
%From  . . .  SLICE  table  1 
Yvh  =  -0 . 005948; 

Nvh  =  -0 . 012795; 

Yrh  =  0  .  001811; 

Nrh  =  -0 . 001597; 

Yvdoth  =  -0.013278; 

Nrdoth  =  -0.000676; 

Yrdoth  =  0  .  000060; 

Nvdoth  =  0.000202; 

%data  extrapolated  from  table  3  of  stability  and  control  derivatives 
for  DSRV  —  Nondimensional 
mfactor  =  -0.97675; 

Izfactor  =  -0.79532; 

m  =  Yvdoth/mf actor ; 

Iz  =  Nrdoth/Izf actor; 

2.  Fin  Specifications  (rudder spec.m) 

%LT  Scott  D .  Roberts 

%Rudder  Specif iciations  for  Towed  Body  USV  retrieval  system 
%Reference  Principles  of  Naval  Architecture  Volume  III 
%20  January  2005 

ar  =  a/2;  %geometric  aspect  ratio 

sweep  =  0;  %sweep  angle  of  the  quarter-chord  line 

Cdc  =  0.8;  %crossflow  drag  coefficient  dependent  on  both  tip 

shape  and  taper  ratio 

Cdo  =  0.0065;  %for  NACA  0015 

e  =  0.90;  %Oswald  efficiency  factor 

taper_ratio  =  0.45; 

bbar  =  croot*ar;  %span  width  in  feet 

ctip  =  croot  *  taper_ratio;  %length  in  feet 
cbar  =  (croot  +  ctip)/2;  %mean  chord  length 

Arud  =  (croot+ctip) *bbar/2 ;  %rudder  area 
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3. 


Total  Hydrodynamic  Derivatives  for  Towed  Body  (hydroderiv.m) 

%LT  Scott  D .  Roberts 

%Total  Hydrodynamic  Derivatives  for  Towed  Body  and  Rudder 
%Reference  Principles  of  Naval  Architecture  Volume  III,  page  242 
%20  January  2005 


%%%%%%%%%%%%  Access  data  for  towed  body  and  rudder 

bodyspec 

rudderspec 

Q-Q-Q-Q-Q-Q-Q.S.Q-Q-Q-Q- 

'o'o'o'o'o'o'o'o'o'o'o'o 


b  =  bbar; 

Alf  =  Arud/ (L*T)  ; 
xFl  =  xF/L; 
bl  =  b/L; 
crootl  =  croot/L; 


%rudder  span  length  in  feet 
%nondimensional  Arud 
%nondimensional  xF 
%nondimensional  b 
%nondimensional  croot 


%Hydrodynamic  Derivatives  for  Rudder  —  Nondimensional 
Yvf  =  -abs (Alf *pi/2*a) ; 

Nvf  =  Yvf*xFl; 

Yrf  =  xFl*Yvf; 

Nrf  =  (xFl ) A2*Yvf ; 

Yvdotf  =  - (4*pi*bl*Alf ) / (aA2  +  l) A (1/2)  ; 

Nrdotf  =  xFl* Yvdotf; 

Yrdotf  =  xFl* Yvdotf; 

Nvdotf  =  (xFl ) A2* Yvdotf ; 


%Total  Towed  Body  Derivates  --  Nondimensional 
Yv  =  Yvh  +  n*Yvf; 

Nv  =  Nvh  +  n*Nvf; 

Yr  =  Yrh  +  n*Yrf; 

Nr  =  Nrh  +  n*Nrf; 

Yvdot  =  Yvdoth  +  n*Yvdotf; 

Nrdot  =  Nrdoth  +  n*Nrdotf; 

Yrdot  =  Yrdoth  +  n*Yrdotf; 

Nvdot  =  Nvdoth  +  n*Nvdotf; 

4.  Stability  Criteria  Calculations  (stability .m) 

%LT  Scott  D.  Roberts 

%Towing  Stability  Criteria  for  Towed  Body 
%Reference :  Latorre  and  Bernitsas 

%28  January  2005 


%Call  data 
hydroderiv 
ro  =  64.1; 
uO  =  8 . 439; 
U  =  uO; 

U1  =  U/uO; 
xp  =  2.4; 


from  hydrodynamic  coefficient  m.file 

%standard  seawater  density  in  lb/ft A3 
%tow  vessel  speed  in  ft/sec 
%towed  body  speed 

%nondimensionalized  towed  body  speed 

%distance  in  feet  from  center  of  body  to  tow  point  at  front 


of  body 
xp 1  =  xp / L ; 

I  =  328; 

II  =  1/L; 

T  =  1000; 


%nondimensionalized 

%length  in  feet  of  unstrained  towline 

%nondimensionalized 

%Tension  in  towline  in  lbf 
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T1  =  T/ ( (ro/2) *LA2*uOA2) ; 


%nondimensionalized 


%Characteristic  equation  coefficients 
A  =  (Yvdot  -  m) * (Nrdot  -  Iz)  -  Nvdot*Yrdot; 

B  =  (Yvdot  -  m) * (Nr  -  Nvdot*Ul)  +  Yv* (Nrdot  -  Iz)  +  (Yvdot*Ul-Yr ) *Nvdot 
-  Yrdot *Nv; 

AO  =  A; 

BO  =  B; 

CO  =  Yv*Nr  +  Nv* (m-Yr) ; 

Cl  =  (m-Yvdot ) *xpl  +  Nvdot; 

C2  =  -(Nrdot  -  Iz)  +  (m-Yvdot ) * (xpl ) A2  +  xpl* (Nvdot +Yr dot) ; 

C  =  CO  +  Tl* (Cl  +  C2/11) ; 

D1  =  Nv  -  Yv*xpl; 

D2  =  -Yv* (xpl ) A2  +  Nv*xpl  +  (Yr+Yvdot) *xpl  -  Nr  +  Nvdot; 

D  =  Tl*  (D1+D2/11) ; 

E2  =  Dl; 

E  =  (Tl /II ) *E2 ; 

alphal  =  (B0*C1*D1  -  A0*(D1)A2)  +  (1/11) * (B0*C1*D2  +  B0*C2*D1  - 
2*A0*D1*D2)  +  (1/ (11) A2 ) * (B0*C2*D2-A0* (D2) A2) ; 
alpha2  =  B0*C0*D1  +  ( 1 /II ) * (B0*C0*D2  -  (B0)A2*E2); 

Gh  =  1- ( (Nv* (Yr-m) ) / (Yv* (Nr-m*xG) ) ) ; 


%Active  Criteria 
R1  =  [xpl  (Nv/Yv) ] ; 

R2  =  [Tl  ( -alpha2 /alphal )] ; 

if  R1 (1)  >  R1  (2) 
condl  =  1; 

else 

condl  =  0; 

end 


if  R2  (1)  >  R2  (2) 
cond2  =  1; 

else 

cond2  =  0; 

end 


5.  Variable  Range  Assignment  and  Data  Collection  (variables.m) 

%LT  Scott  Roberts 

%Program  Run  sequence  for  all  values  of  n,  a ,  croot,  and  xF 
%01  February  2005 

clc, clear 
datav  =  [ ] ; 

datanf  =  [ ] ; 
for  n  =  1  .  0  : 1 : 2 . 0 

for  a  =  0.5:0. 5:6.0 

for  xF  =  -0 . 2 : -0 . 1 : -2 . 5 

for  croot  =  0.1:. 05: 1.5 
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Gh]; 


end 


stability; 

datav  =  [datav;  n  a  xFl  crootl  Rl(l)  Rl(2)  R2(2) 


end 

end 

end 
n  =  0; 

for  croot  =  0.1:0.05:1.5 
stability; 

datanf  =  [datanf;  n  a  xFl  crootl  Rl(l)  Rl(2)  R2(2)  condl  Gh] 

end 


condl 
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